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Transonicity in a spherically symmetric accreting system has been considered in both the sta- 
tionary and the dynamic regimes. The stationary flow, set up as a dynamical system, has been 
shown to be greatly unstable to even the minutest possible deviation in the boundary condition 
for transonicity. With the help of a simple analytical model, and some numerical modelling, it has 
then been argued that the flow indeed becomes transonic and stable, when the evolution of the 
flow is followed through time. The time-dependent approach also shows that there is a remarkable 
closeness between an equation of motion for a perturbation in the flow, and the metric of an analog 
acoustic black hole. 
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I. INTRODUCTION 

o 

Accretion processes involve, in very simple terms, the flow dynamics of astrophysical matter under the external 
gravitational influence of a massive astrophysical object, like an ordinary star or a white dwarf or a neutron star 
£ — | The qualifier "external" is to be stressed upon here, to distinguish accretion processes from the self-gravity driven 
C^*> . collapse of a fluid system, as in the case of a star itself. The accreting astrophysical matter whose fluid properties 
we are interested in, could be the interstellar medium — as modelled by its spherically symmetric infall on to an 
isolated accretor — or stellar matter, as seen in a binary system, where tidal deformation of a star, leads to matter 
flowing out from it into the potential well of a compact companion In all of these cases, the fluid system is 
satisfactorily described by a momentum balance equation (with gravity as an external force) , the continuity equation 
> ■ and a polytropic equation of state. 

In astrophysics, studies in accretion, in a formal sense, have been carried out for more than half a century now. 
• • Initially, astrophysical problems in the nature of what we understand to be accretion processes at present, were 
. £h ' studied by Hoyle and Lyttleton in the context of the infall of matter on to a star moving through the interstellar 
medium [2J. In their methods, however, Hoyle and Lyttleton neglected the pressure effects, with the argument that 
' any heat generated would be radiated away rapidly, so that the temperature of the infalling gas (and related to that, 
the effects of pressure as well) would remain negligibly low [2|. This was found to be a satisfactory prescription for 
most cases of astrophysical interest, which were then being studied. In 1952, however, in a very important paper (2j, 
Bondi attacked this problem somewhat differently, by taking into account the pressure effects. This work on spherically 
symmetric infall of matter, on to a massive and attracting centre, using formal fluid dynamical equations, has assumed 
a paradigmatic status in accretion studies [l|, Q • 

Bondi studied the problem of spherical accretion in its stationary limit, i.e. by only considering the extreme case of 
negligible dynamical effects. In his own words the mathematical difficulties associated with the problem of studying 
both the pressure effects (ignored by Hoyle and Lyttleton) and the dynamical effects were "insuperable" Q, given 
the computational facilities available to him half a century ago. 

The stationary equations would lead to various classes of solutions. Of these the interesting ones would be those 
which obey the outer boundary condition that at large radii the flow velocity would be highly subsonic, i.e. small 
compared with the speed of sound, with the speed of sound itself approaching a constant "ambient" value at large 
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radial distances. With this outer boundary condition being satisfied, it would be physically meaningful to consider 
only the entire class of solutions, which remain subsonic everywhere, and the exceptional case of the lone transonic 
solution. This exceptional solution, determined uniquely by the value of the density at infinity, crosses the sonic point 
of the flow — a point which is given by the radius where the flow velocity smoothly matches the speed of sound — 
and acquires supersonic values at lesser radial distances. Further, in one sense this transonic solution represents a 
limiting value for the mass infall rate, for a given value of the density at infinity, because any inflow rate higher than 
that given by the transonic solution, will make the flow "bounce" outwards Q. 

In this situation it became worthwhile to inquire into the natural selection of a solution by a spherically symmetric 
accreting system. Regarding this question Bondi himself offered some insights when he observed in his paper that 
since it is with the transonic solution that the lowest energy is associated, we may find a natural system in this 
state 0. Furthermore, Bondi suggested that this result would also be "in agreement with the intuitive idea that, 
since there is nothing to stop the process of accretion, it takes place at the greatest possible rate" i-e. the rate 
given by the transonic solution. 

These arguments by Bondi were essentially based on physical considerations. His suggestion that a linear stability 
analysis might offer a clue was also put to the test 0, S Subjecting the physically relevant stationary inflow 
solutions to a linearised time-dependent perturbation led to the result that all solutions displayed stable behaviour, 
and through this method, no insight may be had as to the natural inclination of the system for any particular solution. 
Garlick settled this point quite categorically by stating that "the plausible assumption that subsonic flows are not 
realised in nature, and that critical flows always develop, cannot be justified by a linear stability analysis, but rather 
by the more fundamental arguments given by Bondi" [5[. 

This then was the understanding that became firmly established — that the transonic solution enjoyed primacy 
over the subsonic ones. It was easy to recognise that for black holes, this would be certainly true. For one thing, 
instead of a physical surface, black holes have what is known as an event horizon, which precludes all possibility of 
a pressure build-up at smaller radii; for another, all matter reaching the event horizon must do so supersonically, 
implying that it must display transonic behaviour 

The situation, however, was not so clear-cut if the accretor had a hard surface like a neutron star or a white dwarf. 
For such an accretor, it may be supposed that the accumulated matter would build up pressure near the surface and 
cause the supersonic flow to be shocked down to subsonic levels, although for a neutron star in particular, all accreted 
matter is expected to be efficiently "vacuum cleaned" away, making it easier for the flow to remain supersonic 

Then again there was the conceptual difficulty in understanding the realisability of the transonic solution, by 
only an infinitely precise determination of the outer boundary condition for the stationary equations of the flow. An 
infinitesimal error in the determination of the boundary condition would generate a solution far away from transonicity. 
Within the framework of the stationary picture, it is difficult to imagine that a natural physical system would be so 
precisely tuned. On the other hand, quite intriguingly, this difficulty is resolved when, instead of merely the stationary 
equations, the transonic solution is tried to be numerically generated through a temporal evolution (accounting for 
explicit time-dependence) of the spherically symmetric accreting system [f| G3 • 

In this article we make some attempt in addressing these issues. We take up the stationary flow equations first 
and model them along the lines of the equations governing a dynamical system. This establishes the nature of the 
sonic point of the flow, as that of a saddle point. Quite apart from the fact that a saddle point is inherently unstable, 
among other adverse implications, this analysis illustrates why in the stationary picture, after having started with an 
outer boundary condition, so much difficulty is encountered in generating a solution that passes through the sonic 
point [ll| . The stationary transonic solution is notoriously unstable under even the smallest of deviations from the 
infinitely precise boundary condition that would be needed to generate the solution. 

Since in a realistic situation, an accreting system would have evolved through only a finite span of time, we then 
consider the case of explicit time-evolution of the flow solution(s), starting from a reasonably physical initial condition 
(given at t = 0). This gives us to understand quite convincingly that not only would the transonic solution be selected, 
but that the physical mechanism that is very likely effective in making the selection is the one about which Bondi 
has conjectured (entirely within the confines of the stationary framework) — that the flow solution with the lowest 
specific energy associated with it, would be preferred to all the others. It is important to note that this selection 
mechanism is actually at work through the temporal evolution of the flow, and that it is entirely non-perturbative in 
character [ill ]. 

In keeping with this time-dependent approach to uphold the exclusive selection of the transonic solution from a 
host of various solutions, we have finally carried out an interesting exercise, that is actually perturbative in nature. 
On imposing a linearised time-dependent perturbation on the constant matter flow rate, we have been able to identify 
a remarkable closeness between an equation of motion for the perturbation, and the metric of an acoustic black hole. 
With the aid of this analogy we have argued that contrary to common belief, even a perturbative treatment might be 
conveying a subtle hint in favour of transonicity. 
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II. THE EQUATIONS OF THE SPHERICALLY SYMMETRIC FLOW 

The variables that we need to consider are the velocity (radial velocity only for the spherically symmetric flow), v, 
and the density, p. We ignore viscosity and write down the inviscid Euler equation for v as 

dv dv IdP dV(r) 
Ot Or p Or Or 

where P is the local pressure and V is the potential due to the gravity of the central accretor of mass M, given by 
V(r) = —GM/r. The pressure is related to the local density through a polytropic equation of state P = Kp 1 , in 
which K is a constant, and 7 is the polytropic exponent [13], whose admissible range is given by 1 < 7 < 5/3, with 
this range having been restricted by the isothermal and the adiabatic limits, respectively. To know how p evolves, we 
need the equation of continuity, 

Our system is specified by Eqs.jJTJ) and ©. We are interested in static solutions, the problem of which has been 
defined in Bondi's own words as follows — "A star of mass M is at rest in an infinite cloud of gas, which at infinity is 
also at rest. ..The motion of the gas is spherically symmetrical and steady, the increase in the mass of the star being 
ignored so that the field of force is unchanging" [3]. Since transonic flows are our concern, we require that the static 
flow evolves from v — ► as r — ► 00 (the outer boundary condition) to v > c s (r) for small r, where c s (r) is the speed 
of sound given by (? s = dP/dp = ^Kp 1 ^ 1 . 

The stationary solution implies dv/dt = dp/dt = 0, and hence we have p = p(r) and v = v(r). This requirement 
renders Eqs.([T]) and as 

1, 

dr p dr 



dv IdP GM 
v^- + --J- + -^=0 (3) 



and 

1 dv 1 dp 2 n 

-3- + -/ + -=°; ( 4 ) 
v dr p dr r 

respectively. It is to be noted here that the Eqs.© and ^ remain invariant under the transformation v — ► — v, i.e. 
the mathematical problem for inflows (v < 0) and outflows (v > 0) is the same [l3[ in the stationary state. 

It is in principle possible to eliminate either v or p and solve for the other variable as a function of r. However, 
adopting a slightly different approach, it is possible to recast Eqs.® and in a combined form as 

dr r \ v z — ci I 




FIG. 1: Stationary solutions for spherically symmetric accretion onto a star. The bold solid curves, A and W, represent 
"accretion" and "wind", respectively. The fixed point is at r = ro and (v 2 /c? s ) = 1. A linear stability analysis indicates that 
the fixed point of the flow is a saddle point. The direction of the arrows along the curve A, demonstrates that the transonic 
flow is not physically realisable in the stationary framework. 
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in which we use the speed of sound, c s , to scale the flow velocity. 

The various classes of solutions of Eq.© are shown in Fig. [T] (to be presently read without the arrows). The two dark 
solid curves labelled A and W refer to the accretion flow and the wind flow, respectively. The meaning of "accretion" 
and "wind" flows, which emerges from Fig. |T| is quite clear. At the point v 2 — c 2 , these solutions correspond to a 
finite slope, i.e. a finite value of dv j dr. This implies that when v 2 = c 2 , we must also have 2c 2 = GM/r, implying 
that the intersection point is a critical point. The "wind" and "accretion" flows smoothly pass through the point 
v 2 = c 2 = GM/2r, with the sonic length scale having been labelled as r = tq. 

The question which now arises is that of the natural preference of the accreting system for a particular solution from 
among all the various possible classes of flows shown in Fig. Q] Linear stability analysis of the various stationary inflow 
solutions in real time, indicates that under the influence of a linearised time-dependent perturbation, all solutions 
appear to be stable, and thus it offers no clue as to the selection of any particular solution 0, d, @|- To resolve this 
issue, it would then be very much worthwhile to recall Bondi's conjecture on this point, that the selection would be in 
favour of that solution, with which is associated the least total energy. The transonic branch satisfies that criterion, 
and hence is the choice. 



III. THE SPHERICALLY SYMMETRIC FLOW AS A DYNAMICAL SYSTEM 



Thus far we have obtained a seemingly wholesome picture regarding the realisability of the transonic flow. However, 
upon a closer inspection of Fig. [TJ we find that there is a problem with it. For solutions which represent flows, the 
associated sense of direction has been assigned by an arrow to each solution. An integration of dv/dr would proceed if 
we start with an initial condition v = Uin at r = r; n far away from the star. For a physically realisable flow, an initial 
condition infinitesimally close to a point on the accretion line A, would trace out a curve infinitesimally close to A 
and in the limit would correctly reproduce A, evolving along it and passing through the critical point (sonic point) 
as we integrate dv/dr, obtained from Euler's equation. We will soon show that the arrows on the integration route 
are as shown in Fig. [T] It is obvious from the direction of the arrows here that the stationary spherically symmetric 
transonic accretion flow is not physically realisable, and closely related to this, the fixed point is also seen to be an 
unstable saddle point. 

To obtain any idea about the direction associated with a solution in the steady picture, it is a matter of common 
knowledge in the study of dynamical systems [14] that we cannot turn to Eq.f5|) in its present form. Rather, it would 
be instructive to write Eq.© in a parametrised form 



d_ 



> 2 ) 



2v 2 2c 



dr , , 

dr = ^ V 



GM 



(6) 



As t, which is an arbitrary parameter, evolves (r is not time, since we are dealing with a stationary flow), we generate 
the v(r) curves. The particular curves which represent transonic flow are the curves which pass through the fixed 



point at r = ro, v — vq, obtained from Eq.©, such that 



-sQ 



and 2c 2 



GM/r . The subscripted label 



represents physical quantities at the critical point. We now need to analyse the nature of the fixed point (ro, ^o/ c so) 
Writing v 2 = v 2 + Sv 2 and r = ro + Sr, and linearising in Sv 2 and Sr, we find 



d 

dr 
d 



Sv 1 



(Sr) 



2v 2 



GM 

(7-l)^ 2 -(2 7 -3)-^fr 



2 r 



(7) 



Using solutions of the form Sv 2 
Eqs.([7]) are found to be 



exp (At) and Sr ~ exp (At) , the eigenvalues of the stability matrix implied by 



A = ±4v/2(5-3 7 ). 



(8) 



For the admissible range of 7, i.e. 1 < 7 < 5/3, the eigenvalues are real with different signs and the fixed point 
(ro, 1) in the r — (v 2 /c 2 ) space is, therefore, identified as a saddle point. The arrows (characterising a saddle point) 
which are needed to make our understanding of the stationary phase portrait complete [l4[ , are as shown in Fig. [T] 
The curves which we have labelled "accretion" and "wind" in Fig. [T] are in fact now seen to be the separatrices of 
a dynamical system, and one cannot traverse the length of a separatrix upto the critical point in any finite range of 
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r values. Starting from an initial point to the right of the fixed point on the curve labelled A, one would need an 
infinitely large number of steps to reach the fixed point, and will certainly not cross it. 

To explicitly establish this result we consider the two parametrised equations in Sv 2 and Sr, given by the Eqs.([7]), 
and using them, we write 



d (Sv 2 ) _ d (Sv 2 ) /dT 
d{5r) ~ d (Sr) /dr ' 



(9) 



We then integrate Eq.([9]) in Sv 2 and Sr, and fix the integration constant from the critical point condition, Sv 2 = Sr = 0, 
to obtain, 



Sr- = ~ s0 



-2c 



2 (7 - 1) ± y/2 (5 - 3 7 ) Sr. (10) 



ro (7 + 1) 

Using Eq. (fTU|) in the latter of the two relations given by Eqs.([7|), we get, 

^=±4,^2(5-37) Sr. (11) 

We can integrate Eg. pip , for both roots, from an arbitrary initial value of Sr = [Sr]- ln to a point Sr = e, where e is 
very close to the critical point given by Sr = 0. We thus get 



1 r d(Sr) - 1 1 h 



C 2 0X /2(5-3 7 ) J[Sr ]la Sr c 2 0X /2(5-3 7 ) 



[Sr] 



(12) 



from which it is easy for us to see that for e — ► 0, \t\ — ► 00. This implies that the critical point may be reached 
along either of the separatrices, only after |r| has become infinitely large. That is why the spherically symmetric flow 
cannot be realised. 

It might also be noted from Fig. Q] that within the framework of the stationary picture, there is another obstacle 
of a more practical nature, standing in the way of the realisability of the transonic flow. Each of the inflow solutions 
- the subsonic ones, the transonic one and the so called "bouncing solutions" is to be obtained by its own very 
precisely defined boundary condition. From among the infinitude of possibilites, the boundary condition that would 
exactly reproduce the transonic accretion curve would have to be defined with infinite precision. And yet, even if that 
practical difficulty were to be satisfactorily addressed, by dint of the fixed point being a saddle, the transonic inflow 
solution would still not be generated spontaneously. 

To illustrate all these points further we carry out a simple numerical analysis. An integration of Eq.Q with the 
help of the polytropic equation of state, and its relation to the speed of sound, will give 

- + nc 2 =E, 13 

2 r 

in which, n is the the polytropic index, given by n = (7 — 1) _1 [Hj], while the integration constant is fixed as 
E = nc 2 (00), for the boundary condition v — ► 0, c s — ► c s (oo) for r — ► 00. 
Integration of the continuity equation, as given by Eq.((3]), will yield 

Airpvr 2 — rn, (14) 

in which the integration constant rh (mass accretion rate) is given by [l| 

/ 9 \ (5-3 7 )/2( 7 -l) 

m = vrG 2 M 2 ^_ f \ . (15) 



cf(oo) \5 — 37, 

Combining Eqs. (|13p and (114|) , along with substituting p by its dependence on c s , will finally give 

v 2 ( u \ 1/n GM . 

T +n Ht nc 2 s (oo) = 16 

2 \vr z J r 

in which fi = (m/47rp 00 )c s 2n (cx)). 

We solve Eq. lfTl)]) for v numerically by the bisection method, using the values M — M Q , c s (oo) = 10 km s _1 , 
Poo = 10 -21 kg m~ 3 and n = 2.5. All these values are typical of accretion of the interstellar medium on to an 
average star. Corresponding to a given boundary condition, each value of r in Eq. (|16p would give a set of two real 
and meaningful solutions. The sequence of data points obtained would indicate that the twin solutions would look 
as shown in Fig. [2] (in which the flow velocity has been scaled as the Mach number), which supports our contention 
about the non-realisability of the critical solutions within the framework set up by the stationary equations alone. 
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Radial Distance 



FIG. 2: A plot of the Mach number, v/c s , against the radial distance r (given on a logarithmic scale). The sequence in which 
the data points have been numerically generated has been shown here by the two curves. 



IV. DYNAMIC SELECTION OF SEPARATRICES : A MODEL 



In the Section IIIII we discussed the non-realisability of the stationary transonic solutions if the hydrodynamic 
accretion problem were to be studied solely in the stationary limit. The non-realisability arises due to the fact that 
the stationary transonic solutions are actually separatrices of various classes of solutions (as the direction of the 
arrows in Fig. [T] would show) in a dynamical system. However, it is widely maintained that transonic solutions are 
indeed to be found in a natural system, as has been seen for the case of the solar wind [IH [l6j . which may be 
treated as a spherically symmetric transonic outflow solution. To reconcile this observational fact with our study of 
the stationary flow picture, it must be appreciated that a real astrophysical problem is not stationary but dynamic 
(time-evolutionary) in nature. Therefore we need to take into account explicit time-dependence of the flow variables 
concerned. In that case it becomes evident that the actual velocity profile would not only depend on the radial 
distance, r, but also on time, t. In this dynamic situation, it is our contention that all the adverse implications 
regarding transonicity would disappear, and a transonic flow would be realised. 

To assure ourselves that such indeed should be the case, we first consider a tractable mathematical model problem. 
We choose a differential equation given by 



dy_ 

dx 



whose integral can be written as 



y 



Ax + 2xy = -C, 



(17) 



(18) 



with C being a constant. This is the equation of a hyperbola. If we want that particular solution which passes through 
the point where f(x, y) = g(x, y) = 0, namely x = y = 1, then C = 2. The curve x 2 — y 2 ~ 4x + 2xy — —2 factorises 
into a pair of straight lines : y — x(l + \/2) + \/2 = and y — x(l — \/2) — \[2 = 0, which are the asymptotes of the 
hyperbola. This pair is shown in Fig. [3] as the lines marked A' and W. 

We want to explore the process of drawing the line A' from a given starting condition. On the line A', we have 
y = at x = 2 + \/2. Let us begin our generating a solution with the condition y = 0atx = 2 + \/2 — e, where 
< e € 1. This starting condition fixes the constant C as C = 2[1 + \/2e — e 2 /2]. Using this value of C, we can 



plot the curve given by Eq. (|18[) . For a given value of x, the value of y is given by the relevant root of the quadratic 
equation thus obtained. The two roots are 



V 



= x ± V2 \(x - if + V2e - e 2 /2 



1/2 



Clearly, to satisfy y — 0atx = 2 + \[2 — e, the negative sign has to be chosen in Eq. (TlT)]) . which will give 

1/2 



y = x 



V2 (x - l) 2 + V2e - e 2 /2 



(19) 



(20) 
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FIG. 3: Integration of Eq.JTTJ) gives a pair of straight lines, with the integration constant fixed by the intersection point (1, 1). 
In the figure the lines are marked A' and W'. Linear stability analysis indicates that for A' and W', the intersection point 
(1, 1) is actually a saddle point, for which the arrows are as shown above. 



At x = 0, y = -V2(l + V2e- e 2 /2) 1/2 , very different from y = which one gets on the line A'. In the limit of 
e — > 0, one generates a part of A' (x > 1) and a part of W (x < 1), instead of the whole of line A'. Another way of 
stating this is that the tracing of A' is utmostly sensitive to initial conditions. If we make an error of an infinitesimal 
amount e in prescribing the initial condition on A', i.e. if we prescribe y = Qa,tx = 2 + \[2 — e instead of y = at 
x = 2 + \/2, then the "error" made at x = relative to A' is 2^/2 which is 0(1). An infinitesimal separation at one 
point leads to a finite separation at a point a short distance away. This is what we mean by saying that the line A' 
(and similarly W) should not be physically realised. 

The clearest and most direct understanding of the difficulty is achieved by recasting Eq. (fl~7|) as a first-order au- 
tonomous dynamical system, described by the set of differential equations 

dy . „ 

dr = x + y ~ 2 

dx .„ . 

- = „-* (21) 

with t being some convenient parametrisation. The fixed point of this dynamical system is (1, 1), namely the point 
where f(x, y) and g(x, y) vanish simultaneously — the point through which A' and W pass. Linear stability analysis 
of this fixed point in the r parameter space, shows that it is a saddle point, with the eigenvalues A given by A = ±V2. 
The solutions passing through the critical point in this x — y space can now be drawn with arrows and the result is 
as shown in Fig. [3J The distribution of the arrows, characterising a saddle point [L4|, implies A' and W cannot be 
physically realised. 

We now investigate if these apparently non-realisable separatrices in the stationary limit may indeed be realised 
when we follow the evolutionary dynamics of y through another variable t. Accordingly, we return to our pedagogic 
example of Eq. (|17|) but now consider y as a field y(x,t) with the evolution through t represented by 

% + { y- x) °l=y + X -2. (22) 

The stationary solution y(x) satisfies Eq.([T7]) and the discussion that follows Eq. (fTT|) is valid for y(x) here. The 
stationary solutions y(x) are as shown in Fig. [3] and the separatrices are y(x) — x(l + \/2) — \[2 and y(x) — 
—x{s/2 — 1) + \/2. We will now show that the dynamics actually preferentially selects these separatrices. 

The general solution of Eq. ([22|) can be obtained by the method of characteristics [13] ■ The two characteristic 
solutions of Eq. (|2"2"|) are obtained from 



dt dx dy 
1 y — x y + x — 2 

They are 



(23) 



y 2 - 2xy - x 2 + Ax = C 



1 , 

x - 1 ± -= {y - x) 



C (24) 
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with the latter having been derived by integrating the dx/dt equation with the help of the solution of the dy/dx 
equation. The general solution of these two equations is given by the condition C — C(C); where C is an arbitrary 
function, whose behaviour is to be determined by the initial condition at t = 0. 

As in a physically realistic situation, we impose the condition that the evolution is driven through a positive range 
of values of t ("time"). This requirement will choose the upper sign in Eqs. ([24| . and the general solution of Eq. ([22|) 
can then be written as 



y 



2xy - x- + Ax = ( 



1 + T2 iy ' x) 



We give the initial condition that y(x) = at t = for all x. This leads to 

1 



C 



x[l- 



V2 



- 1 



4x, 



and gives a form for the function £ as £(z) = Az 2 + Bz + C, in which, 



-V2t 



(25) 



(26) 



A= — 



B = 



With the initial condition y(x) = at t = 0, the solution to Eg. ([23)) reads 



C = 2. 



y-x 



V2 + 1) + V2 



y + x 



V2- 1) - V2 



with 



,y) = x[i 



1 

71 



i 



B4>e 



y 

V2~ 



-V2t 



Atfe- 2 ^ 



(27) 



Clearly as t — > oo, the right hand side in Eq. ([27[) tends to zero and we approach one of the two separatrices (which 
were otherwise non-realisable from the stationary viewpoint) shown in Fig. Of the two separatrices, the one which 
will be relevant will be determined by some other imposed requirement. For the astrophysical flow, the two separatrices 
are the transonic accretion and the wind solutions. One chooses the proper sign of the velocity (v < for inflows, 
and v > for outflows) to get the flow in which one is interested. 



V. A NON-PERTURBATIVE SELECTION OF THE TRANSONIC FLOW 



We have seen with the help of a model problem that the dynamics makes it possible as a matter of a mathematical 
principle, to select solutions which were apparently non-realisable from the stationary perspective. We now extend 
that treatment to the accretion problem here. It is important to understand that a real astrophysical flow itself is 
dynamic in character. This implies that explicit time-dependence of the flow equations would have to be taken into 
account. Having said that, it must also be said that the Eqs.{TJ) and ([2]), which govern the temporal evolution of the 
flow, do not lend themselves easily to a ready mathematical analysis; indeed, in the matter of incorporating both the 
dynamical and the pressure effects in the equations, short of a direct numerical treatment, the mathematical problem, 
in Bondi's own word — "insuperable" 0] — is very appropriately described. Therefore, to have an appreciation of the 
governing mechanism that underlies any possible selection of a transonic flow, we have to adopt some simplifications. 

In our accretion problem we study the dynamics in the regime of what is understood to be the pressureless motion 
of a fluid in a gravitational field [18| — which is a line of attack that is somewhat reminiscent of the methods of 
Hoyle and Lyttleton, as Bondi has mentioned in his paper Simplification of the mathematical equations, however, 
is not the only justification for such a prescription. A greater justification lies in the fact that the result delivered 
is in conformity with, what Garlick calls "the more fundamental arguments of Bondi" Q, that it is the criterion of 
minimum total energy associated with a solution, that will accord it a primacy over all the others. 

An immediate consequence of adopting dynamic equations is that the invariance of the stationary solutions under 
the transformation v — > —v, is lost. As a result, we now have to separately consider either the inflows (v < 0) or 
the outflows (v > 0), a choice that we impose upon the system at t — 0. Eulcr's equation, tailored according to our 
simplifying requirements, is rendered as 

dv dv GM „ 

m +V 8-r + — =°' (28) 
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which we solve by the method of characteristics [17[ • The characteristic curves are obtained from 

dv 



dt 

T 



dr 

v 



-GM/r 



On first solving the dv / dr equation, we get 



v 
~2 



GM 



2 ' 



(29) 



(30) 



with c being an integration constant obtained from the spatial part of the characteristic equation. We use this result 
to solve the dr/dt equation from Eq.([29|). and for c 2 > 0, we get 



- eft) - In 



(31) 



in which c is another integration constant, and r s is a length scale in the system defined as r s — 2GM/c 2 . A similar 
expression can also be written for c 2 < 0. 

A general solution of Eq. (|29p is given by the condition, c = £(c 2 /2), with £ being an arbitrary function, whose form 
is to be determined from the initial condition. We can, therefore, set down the general solution as 



— (yr — c 2 i) — In 



r / v \' 
- - + 1 

7% V C / 



v 2 GM 



(32) 



to determine whose particular form we use the initial condition, v — u^{r) at t — for all r, where uq is in general 
some initial velocity distribution over space. It should be easy to see that for t — > oo, we would get the stationary 
solution 



GM 



= 0, 



(33) 



with the long-time evolutionary approach towards this stationary state behaving as i~ 2 / 3 . 

For a simple intuitive understanding of the physical criterion that drives the flow towards a chosen stationary end, 
we set Mo = 0. This initial condition will necessitate c 2 < 0. The whole physical picture could be conceived of as one 
in which a system with a uniform velocity distribution v = everywhere, suddenly has a gravity mechanism switched 




FIG. 4: Evolution of the velocity field as given by Eq. (|29[) , under the initial condition v = at t = for all r. The horizontal 
line on top of the plot represents a limiting value for the velocity, \J2GM jr , which is being terminally approached by —v, 
whose evolution through t is being followed at the fixed length scale, r — 51vq, with tq being the radius of the accretor and 
with M = M Q . 
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FIG. 5: The slope of this logarithmic plot shows that in the early stages of the evolution, —v varies linearly with t. Deviation 
from this linear growth sets in on time scales of 10 5 seconds. The horizontal line on top is the limiting value of the velocity. 



on in its midst at t — 0. This induces a potential —GM/r at all points in space. The system then starts evolving to 
restore itself to another stationary state, so that for t — > oo, the total energy at all points, E = (v 2 /2) — (GM/r) = 0, 
remains the same as at t = 0. This is evidently the stationary solution associated with the lowest possible total energy, 
and the temporal evolution selects this particular solution from all possible meaningful solutions. 

This result has been borne out by a numerical integration of Eq. (|29[) by the finite differencing technique. The mass 
of the accretor has been chosen to be M© , while its radius is r Q . The evolution through time has been followed at 
a fixed length scale of 51r©. The result of the numerical evolution of the velocity field, — v (for inflows v is actually 
negative), through time, t, has been plotted in Fig. [5] The limiting value of the velocity, as the evolution progresses 
towards the long-time limit, is evidently y/2GM/r (with M — M Q and r = 51r©), as Eq. ([33|) would give us to believe. 
This is what the plot in Fig. 2] shows, as — v approaches its terminal value for t — > oo. The slope of the logarithmic 
plot of —v against t in Fig. [5] indicates that in the early stages of the evolution there is a linear growth of the velocity 
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FIG. 6: Evolution of the velocity field, scaled as the Mach number, M, through time, t. Transonicity is clearly evident, as all 
curves cross the M = 1 line. Moving from left to right, successive solutions have been shown for t = 1000, 2000, 3000 and 4000 
seconds, respectively. The radial distance along the horizontal axis has been scaled by the sonic radius, ro. 
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FIG. 7: Progression of the sonic front in the early stages of the temporal evolution of the flow. Once again the radial distance 
along the horizontal axis has been scaled by the sonic radius, ro. The slope indicates a power dependence of 1/2 for r versus t. 



field through time, but on later times, conspicuous deviation from linearity sets in. 

The foregoing argument can now be extended to understand the dynamic selection of the transonic solution. The 
inclusion of the pressure term in the dynamic equation, fixes the total energy of the system accordingly at t = 0. A 
physically realistic initial condition should be that v — at t — 0, for all r, while p has some uniform value. The 
temporal evolution of the accreting system would then non-perturbatively select the transonic trajectory, as it is this 
solution with which is associated the least possible energy configuration. This argument is in conformity with Bondi's 
assertion that it is the criterion of minimum total energy that should make a particular solution (the transonic solution 
in this case), preferred to all the others. However, this selection mechanism is effective only through the temporal 
evolution of the flow. 

To test this contention a numerical study has been carried out, once again using finite differencing, but this time 
using both the dynamic equations for the velocity and the density fields, as given by Eqs.([T]) and ([2|). The accretor 
has been chosen to have a mass, M , and radius, r Q . The "ambient" conditions are c s (oo) = 10 km s^ 1 and 
Poo — 10~ 21 kg m~ 3 , while the polytropic index, n = 1.6. For these values of the physical constants, transonicity 
becomes apparent even at the very early stages of the evolution. The course of the evolution of the velocity field 
(scaled by the speed of sound), at various points of time, for a substantially representative range of the radial distance 
(scaled by the sonic radius) has been shown in Fig. [5] 

The outward propagation of the sonic front, as time progresses, has been traced in Fig. [7j In this logarithmic 
plot, what can be seen in the early stages of the evolution, is that the sonic front travels through space with a 1/2 
power dependence on time. This behaviour can be set down as r = roQi 1 / 2 , with the cons tant factor Q having been 
determined empirically from dimensional considerations as Q — 4(5 — 37) _1 \J c^(oo) / G M . This estimate tallies very 
closely with the numerical value obtained from the plot in Fig. [Jj nothwithstanding which, a cautionary note that has 
to be sounded here is that this quantitative match holds good for the early stages of the evolution only and need not 
be strictly applicable for the entire span of the temporal evolution of the velocity field. 

VI. A PERTURBATION EQUATION AND THE METRIC OF AN ACOUSTIC BLACK HOLE 

We have discussed many times in the previous sections that subjecting the stationary solutions to a linearised 
time-dependent perturbation shows that all the acceptable solutions are stable. Therefore, through a perturbative 
analysis, no clue could be had as to the special status of any solution. We subject this line of thinking to a closer 
inspection. 

To carry out a linear stability analysis Petterson et al. [J] and Theuns and David [6j have made use of a variable 
defined as / = pvr 2 , whose stationary value is to be obtained from the continuity equation, given by Eq.Q, and this 
background value, /t>, is seen to be a constant that is physically identified with the mass flux. In spherical symmetry, 
the flow variables are v and p. If we impose small perturbations, v' and p' , on the stationary background solutions, 




12 



Vb and pb, we may derive a linearised relation for the perturbation of / as /' = (v'pb + Vbp') r 2 , and in terms of this 
perturbed quantity we may then obtain a linearised equation of motion for the perturbation, given by 



0_ 

at 



^ (d_r_ 

fb V dt 



d_ 

dt 



d_r_ 

dr 



d_ 

dr 



fb 



df_ 
dt 



d_ 

dr 



Vb_ 

fb 



= 0, 



(34) 



with c s b being the background stationary value of the speed of sound. 

We now approach this whole question from a different perspective. It is known that there is a close one-to-one 
correspondence between certain features of black hole physics and the physics of supersonic acoustic flows [19 '] . For 
an irrotational, inviscid and barotropic fluid flow, Euler's equation may be written as 



dv 1 VP 

^+ Vv.v + — + W = 0, 

dt 2 p 



(35) 



in which V ~ —GM/r, as before. In this situation we can represent velocity as the gradient of a scalar function ip, i.e. 
v = —Vip. For the barotropic condition p = p(P), it is possible to write V/i = (VP)/p, upon which, from Eq. 
we may derive the result 



d ^ , 1 m i\ 2 

-at + 2 (w) 



h+ V = 0. 



(36) 



On a time-dependent background solution (pb, Pb,ipb) we impose a perturbation (p' , P',ip'), and together with Eq.([36 
and the continuity equation given by Eq.Q, a linear equation for the perturbation is delivered as 



d_ 

at 



Pb_ 

r 2 

C sb 



dip' 

~at 



+ v b • Vip' 



+ V 



-PbW 



Pb 

<4 



dip 1 

~at 



+ v b • Vip' 



0. 



(37) 



in which cj h = P 1 j p' . This result represents an equation of motion of an acoustic disturbance. Visser discusses this 
subject [19| by stating that if a "fluid is barotropic and inviscid, and the flow is irrotational (though possibly time- 
dependent) then the equation of motion for the velocity potential describing an acoustic disturbance is identical to the 
d 'Alembertian equation of motion for a minimally coupled massless scalar field propagating in a (3 + 1) - dimensional 
Lorentzian geometry." In other words, Eq. (|37[) may alternatively be represented by a compact formulation given as, 



d 
dxf* 



(38) 



in which the Greek indices run from to 3. The two expression given by Eqs.(|37)l and (.38) arc completely equiva- 
lent [lj| . If we consider the special case of only radial dependence for g^ v , then identifying t as and r as 1 in the 
Greek indices, we will get from Eg. ([38]) . 



.g 00 = l, g^=g w 



Vb, 



g 11 =v 2 - c 2 sb . 



(39) 



Interestingly enough, the expression for /', as Eq. ([34| shows it, can also be reduced to a similarly compact form 
as Eq. (f3"5)) gives for ip', with an equivalent identification for an effective metric g**" . It is easy to compare terms 
between Eas. (|34]) and ([38]) with ip' replaced by /', and see that the same set of g^ v is obtained for /', as Eq. (|39"l) 
shows for ip' . This closeness of form is very intriguing. The physics of supersonic acoustic flows closely corresponds 
to many features of black hole physics. For a black hole, infalling matter crosses the event horizon maximally, i.e. 
at the greatest possible speed. By analogy the same thing may be said of matter crossing the sonic horizon of a 
spherically symmetric fluid flow, falling on to a point sink. That this fact can be appreciated for the spherically 
symmetric accretion problem, through a perturbative result, as given by Eq. p4[) . is quite remarkable. This is because 
conventional wisdom tells us that we would be quite unable to have any understanding of the special status of any 
inflow solution solely through a perturbative technique 0, H| . It is the transonic solution that crosses the sonic horizon 
at the greatest possible rate, and the similarity of the form between Eqs. flM)) and (|3"7]) may very well be indicative 
of the primacy of the transonic solution. If such an insight were truly to be had with the help of the perturbation 
equation, then the perturbative linear stability analysis might not have been carried out in vain after all. 



VII. CONCLUDING REMARKS 



We have seen how achieving transonicity becomes distinctly easy when one accounts explicitly for time-dependence 
in the mathematical problem of spherically symmetric accretion. Certainly the physical and mathematical difficulties 
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associated with a purely static approach to this question, disappear immediately upon involving time. But that is 
not to say that all questions have been answered satisfactorily. For instance, a very important issue that has to 
be addressed in greater detail is the manner in which the temporal evolution drives the velocity field towards its 
stationary critical (i.e. transonic) end, especially in the long-time limit. Even in the simplified pressure-free regime, 
we have seen for ourselves that this is not something whose answer may be given very easily. With the involvement 
of the evolution of both the velocity and the density fields — as it has to be for a real astrophysical problem - 
the computational difficulties will be quite staggering. Nevertheless, this has to be the subject of a more intensive 
scrutiny. 

A further intriguing issue is whether or not a perturbative analysis in real time — considered of not much help in 
understanding the primacy of the transonic state — can indeed offer some insight into questions related to transonicity. 
The closeness of an equation of motion for a perturbation in the accretion problem, to the metric of an acoustic black 
hole, has been a beguiling revelation in this regard. 
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